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Abstract 

The inflaton must convert its energy into radiation after inflation, which, in a 
conventional scenario, is caused by the perturbative inflaton decay. This reheating 
process would be much more complicated in some cases: the decay products ob- 
tain masses from an oscillating inflaton and thermal environment, and hence the 
conventional reheating scenario can be modified. We study in detail processes of 
particle production from the inflaton, their subsequent thermalization and evolution 
of inflaton/plasma system by taking dissipation of the inflaton in a hot plasma into 
account. It is shown that the reheating temperature is significantly affected by these 
effects. 



1 Introduction 



The idea of inflation [H [2] has now become a part of standard cosmological evolution 
scenario. It provides beautiful explanations for the nearly flat isotropic/homogeneous 
Universe and the origin of primordial density fluctuation, which results in rich observed 
cosmological structures. For successful inflation, the energy of the inflaton, which drives 
the inflationary expansion of the Universe, must be transferred to the radiation consisting 
of hot standard model (SM) plasma. This process, called reheating, is rather an unknown 
aspect of inflation [5], partly because the process is model dependent and partly because 
the era of reheating is difficult to be explored observationally. However, the reheating 
temperature Tr, corresponding to the temperature of the hot plasma at the beginning of 
the radiation dominated era, is an important characteristics of inflation model since it often 
determines the efficiency of leptogenesis/baryogenesis and the abundance of (unwanted) 
relics such as the gravitino and moduli. 

In a conventional picture, the inflaton is assumed to have a coupling to light fields 
and perturbatively decays into them. Produced light SM particles are thermalized and 
constitute radiation component of the Universe thereafter. In this case, the reheating 
temperature simply depends on the perturbative decay rate of the inflaton: 



where is the relativistic degrees of freedom at the temperature T = T^' ' , r° denotes 
the inflaton decay rate evaluated at the vacuum and Mp\ is the reduced Planck scale. Here 
we defined this quantity, T^'^ \ as "would-be-reheating temperature". 

However, this simple picture does not hold for some inflation models for the following 
reasons. First, the inflaton is oscillating around its potential minimum, and hence masses 
of coupled particles also oscillate with time, which would invalidate the use of the inflaton 
decay rate at the vacuum [H El [6]. Second, before the complete decay of the inflaton, 
the Universe is often already filled with high-temperature plasma. Thus light particles, 
including SM particles, obtain thermal masses and should be treated as quasi-particles, 
which would significantly modify the inflaton decay rate into these particles. On this 
second aspect, it was pointed out that large thermal masses of decay products prevent 
the inflaton decay and the temperature of plasma cannot be as high as the inflaton mass 
(divided by a coupling constant) [7]. This is not true: in high-temperature environment, 
the quasi-particles obtain thermal widths and the inflaton dissipates its energy into thermal 
plasma as was explicitly shown in Refs. [8^ i9j in the context of reheating. Intuitively this is 
understood as a result of efficient scattering processes between inflaton and quasi-particles 
in thermal plasma. Thus actual thermal history would be much more complicated and, 
in particular, the reheating temperature would be significantly different from the estimate 



In this paper we address the issues of thermalization and reheating after inflation in 
detail. We start from the inflaton oscillation just after inflation and show how particle 
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production and their thermalization occur. Then we study the evolution of the inflaton 
oscillation and plasma by taking into account the inflaton dissipation in high-temperature 
plasma and also the non-perturbative particle production, until the inflaton dissipates all 
its energy, after which the radiation dominated Universe begins. Formulations for these 
effects are found in our previous paper [10], where dynamics of scalar fields in thermal 
environment was studied in detail. 

It is shown that deviation from the conventional reheating scenario becomes more 
prominent for smaller inflaton mass, and hence we mainly focus on low-scale inflation 
model. Such low scale inflation is actually realized in the Higgs inflation ^TTl [12] , where 
the SM Higgs field plays a role of inflaton, since its mass around the vacuum is weak scale. 
Some supersymmetric (SUSY) inflation models are also classified into this category: e.g., 
MSSM inflation [T3], alchemical inflation [H] or others. 

In Sec. 121 we briefly discuss particle production and their thermalization. Dissipation 
coefficients in thermal plasma are also listed. Using these ingredients, we study the evo- 
lution of inflaton and plasma system, and determine the reheating temperature in Sec. [21 
We conclude in Sec. [H 

2 Particle production and dissipation 
2.1 Setup 

Let us consider a following simple setup where the inflaton interacts with light fields x 
via Yukawa interaction: 

C = £kin - + A0 (XlXR + h.c.) + £other (2.1) 

where A is a coupling constant taken to be real and positive, £kin denotes canonical kinetic 
terms, and £other denotes the other light degrees of freedom including gauge bosons. The 
bare mass of x is neglected in what follows. We also assume that the x fields are charged 
under some gauge groups and they interact with other light degrees of freedom via these 
gauge interactions. The coupling constant A and gauge coupling g are assumed to be 
smaller than unityl^ We also define a = g"^ /{An) for later convenience. 

Note that the model (12.11) should be regarded as a representative model which correctly 
describes essential features of more general class of models. It is straightforward to extend 
the model as 

^ = >Ckin - Y \^lkA + Yl ^''^^^^ {XL,kXR,l + h.c.) + £other, (2.2) 

It is possible to consider the case where Yukawa interactions doniinantly connect the x fields with the 
other light degrees of fi'eedom. The following calculation does not change much if the coupling constant 
g is reinterpreted as the Yukawa coupling. 



2 



where integers k, I include possible flavor and gauge indices. In particular, the inflaton 
can have gauge charges. If one of the scalar fields (pki obtains a large field value and takes 
a role of inflaton, the dynamics of inflaton is effectively described by a simplified model 
(12. ip . In a SUSY model, there should be a inflaton coupling to bosons as C = A^|0p|xP 
with X denoting the scalar partner of x- Inclusion of this coupling does not modify the 
following arguments, as long as we restrict ourselves to the case where the parametric 
resonant phenomena do not occur (see Sec. 12. 2p . 

The initial amplitude of inflaton at the end of inflation, (pi, is taken as a free parameter, 
without specifying the inflaton potential beyond (pi. We only assume that the subsequent 
coherent oscillation can be well described by the quadratic potential: m^0^/2. As men- 
tioned in the Introduction, for the small amplitude case, i.e. ^ A0 with representing 



the amplitude of oscillating inflaton fieldj^j the effect of thermal plasma on the inflaton 



dissipation was already studied by [H [9] . Hence, hereafter, we mainly concentrate on the 
large amplitude case: -C A0. The purpose of this paper is to clarify the thermalization 
and reheating process in this class of models, especially the reheating temperature, by 
solving the evolution of inflaton and plasma system. 

The basic ingredients that we will use in the following discussion are found in detail in 
our previous paper [TU] and here we briefly repeat the results. Let us see what happens 
after inflation in the following subsections. 



2.2 Particle production and thermalization 
2.2.1 Instant preheating 

After the inflation, the inflaton starts to oscillate around its potential minimum with an 
initial amplitude 0j. Hence, the coupled field x has an amplitude-dependent dispersion 
relation: uj"^ = k'^ + m^^(Ty + A^0^(t) where m^^{T) denotes a possible thermal mass. 
Initially the thermal mass vanishes {m^^ = 0), since there is no background plasma. If 
X(f)i <^ m^, the perturbative decay of the inflaton into x creates thermal plasma as in a 
conventional scenario. On the other hand, if A0j ^ m^, such a process is kinematically 
blocked due to the 0-dependent mass of x at the most time domains in each one oscilla- 
tion of the inflaton. Instead, the following non-perturbative particle production process 
becomes important J^^ 

The efficient non-perturbative particle production occurs when the adiabaticity of the 
coupled X fields is broken down [1]: IcO^/w^l ^ 1. From this inequality, the following 
condition is obtained: 

X(f) ^ max 



(p and (j){t) represent the amplitude and oscillating field value respectively. 
T^"^ The perturbative inflaton decay into gauge bosons through one-loop process involving x held is 
possible, but its efficiency is lower than that of the non-perturbative particle production. 



(2.3) 
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Here stands for the amplitude of 0. Importantly, Eq. (12 .Sp implies that if the thermal 
mass m^jj(T) ~ gT becomes as large as fc* = (A^m^)^/^, then the non-perturbative particle 
production does not occur [10]. If Eq. ( 12. 3p is met, the x's modes below the typical 
momentum are amplified in each oscillation, and the typical number density for one 
degree of freedom can be evaluated as 

^ (2^0^ ~ 8vr3 ' (2.4) 

in each oscillation. 

After the (j) passes through the origin, the produced particles become heavy due to 
the field value of oscillating 0, and their decay rate becomes large correspondingly. If 
the decay rate of x'- sufficiently large, they eventually decay into the other hght 

degrees of freedom at r^(0(tdec)) ec ~ 1 well before the reaches its maximum value [6] . 
This is the case for <^ aX(j), if we assume that the typical decay rate of x is given by 
~ am^ ~ aA|0()f:)|. In this linear potential for generated by the produced 

particles |4j is insignificant, since the produced particles decay when the infiaton field value 
reaches ~ [m^cj)/ {aX)Y^^ <^ 0. 

In addition, if the infiaton couples to the bosonic fields x, this condition guarantees the 
absence of violent parametric resonant phenomena. Unless the coupled bosonic particles 
decay before the infiaton moves back to the origin, the production rate of these bosonic 
particles is enhanced due to the induced emission effect. Thus, their number density 
grows exponentially whenever the infiaton passes through the origin, and as a result, the 
system may enter the so-called turbulent regime as shown in Refs. |T5] with classical lattice 
simulations and recently in Refs. [16] with Kadanoff-Baym eqs. 

Throughout this paper, we assume that the condition <^ aA0 holds, and hence 



such effects can be neglected]^ The energy density converted to the other light degrees 



of freedom in one infiaton oscillation can be evaluated as 

6p ~ m^n^l^^^ ~ a"^/^(A0m0)^ (2.5) 



2.2.2 Thermalization of plasma 

In order to study the subsequent evolution of the produced plasma and oscillating scalar 
field at every moment, it is practically important to know whether or not the produced 
other light degrees of freedom can attain thermal equilibrium in a time scale faster than 
the oscillation time scale of |17] . 

At the first passage of ~ 0, the total energy density of light degrees of freedom prad 
is given by prad = Sp, since there are no particles before this non-perturbative particle 

The number density of particles produced by the decay of x fields directly coupled to the infiaton 
cannot become so large since otherwise the non-perturbative production becomes inactive due to the large 
screening mass: mg > fc* [cf. Eq. (|2.3|) ]. 
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production. In this case, as extensively discussed in Ref. [18], the thermahzation time 
scale of light degrees is estimated by 



teq ~ (a'Tf)"' ^/QjTf ~ «-33/l6(A0m^)-l/^ (2.6) 

where Tf ~ p^J^ and the typical momentum scale Q is given by Q ~ "^xldec ~ ^"^''^(A^m,^)^/^. 
[See Appendix |X] for more detail.] Therefore, the produced other light degrees of freedom 
can be safely regarded as "thermal" plasma as far as the following condition is met: 

l<^a^^'^\Mlm^fl\ (2.7) 

As one can see, this condition is satisfied if the initial amplitude of oscillating scalar field 
is sufficiently large. 

After several oscillations, the energy density of background thermal plasma becomes 
much larger than the one produced via the non-perturbative production in each oscilla- 
tion, i.e. pra.d ^ Sp- In this case, the equilibration time is given by the relaxation one, 
corresponding to a time scale for a hard particle Q > T to emit its energy away to the 
thermal plasma of temperature T 



#5 



trix- («'T)"iv/g7T. (2.8) 
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Hence, the background plasma can remain in thermal equilibrium if trix ^ 

If these conditions are met, the produced light particles attain thermal equilibrium in 
each oscillation, and consequently the screening mass of coupled x field can be described 



by the thermal mass ~ gTi^ As mentioned above [Eq. (12. 3p ]. the instant preheating 
stage finishes when this thermal mass becomes comparable to ~ (A^m^)^/^. 



2.3 Dissipation to thermal plasma 

As discussed in the previous subsection, the background thermal plasma is produced via 
the instant preheating if A0j > m^, just after infiation. The produced thermal plasma 
can significantly affect subsequent dynamics of oscillating inflaton field. Aside from the 
blocking effect on the non-perturbative production due to the thermal mass, which we 
discussed in the previous section, there are basically two effects from thermal environment: 
(i) thermal effective potential for the inflaton and (ii) dissipation of the inflaton to thermal 
plasma. We are mainly interested in the situation where the scalar fleld dominates the 
Universe, and hence let us discuss the latter effect (ii)cl 

Note that this temperature T is not related to the Q, contrary to the Tf. 

Otherwise, the effective mass for x cannot be described by a temperature T and following analyses 
become more complicated. We do not go into such a case in this paper. 

If the scalar field 4> oscillates dominantly with the thermal potential {e.g. thermal mass or thermal 
log), its energy density is bounded as < [10]. Therefore, it is typically less than the energy density 
of thermal plasma: prad ^ g*T^ > P(t>- 
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We will not perform detailed calculations of the dissipation coefficient and will not 
show all the list of dissipation coefficients in various regimes in this section. Instead, let us 
explain its typical behavior and intuitive physical interpretation relevant to our following 
discussion. We refer to Refs. P UHl [iHl EH] for basic formalism to calculate them. The 
complete list of dissipation coefficient in various regimes is summarized in Appendix |Bl 

The effect of thermal plasma becomes significant in the case where the typical time scale 
of oscillation is much slower than that of interaction in thermal plasma, i.e. <^ aT. 
In this case, the oscillating scalar field cannot decay into light degrees of freedom, 
since these would-be decay products acquire thermal masses which are larger than the 
mass. However, the oscillating scalar can dissipate into thermal plasma through multiple 
scattering by light particles in thermal plasma, or more precisely through thermal width 
of each quasi-particle excitation. 

The dissipation coefficient depends on the value of 0, and its dependence can be divided 
into two regimes: (i) small field value regime A0 <^ T and (ii) large field value regime 
A0 ^ T. In the case (i), the coupled x particles are relativistic and its number density 
is given by T^. Hence, the oscillating scalar dissipates its energy through scatterings 
involving x particles. On the other hand, in the case (ii), such processes are unlikely 
to occur since the x particles become very heavy due to the amplitude of and the 
number density of x is exponentially suppressed correspondingly. Therefore, the oscillating 
scalar dissipates its energy mainly by multiple scattering of gauge bosons through a higher 
dimensional operator obtained from integrating out the heavy x field. Consequently, the 
dissipation coefficient for nis ^ oT can be evaluated as [10] 



' Aq dim(r)A2ar/(27r2) for A0 < ~ gT 
Ao dim(r)A V/(vr^aT) for ^ gT ^Xcjx^T (2.9) 
ba^Ty(j)^ for A0 > T 



where 



/T(r)\' (127r)2 ^ ^ 

Here dim(r) stands for the dimension of x's representation r of gauge group and T(r) is 
the index of x's representation r defined by T{r)6°-^ = tT[t°-{r)t^{r)], and Aq is a numerical 
constant, typically Aq ~ 1/2. For our numerical calculation, we take a = 0.05, and then 
it is given by Aq ~ 0.3. Note that the above dissipative coefficient is calculated in two 
limits: large and small amplitude, and hence we have some ambiguities in the intermediate 
regime, 
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In the opposite limit, > T, the dissipation coefficient can be estimated with ne- 
glecting the finite density correction to the dispersion relation of x- Therefore, it is simply 



In addition, the small amplitude result computed with one-loop approximation may change by some 
factors due to the resummation of infinitely many higher- loop diagrams as discussed in Ref. |20] . 
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given by the perturbative decay rate of oscillating scalar into light degrees of freedom. 
If the amplitude is much larger than the mass of {i.e. A0 ^ ^(f>), the oscillating 
scalar loses its energy mainly via the non-perturbative particle production as discussed 
in the previous section. Hence, practically, the pertrurbative decay becomes important at 
A0 -C and it is given by 



dim(r) 



6TT 



{2.11] 



In the intermediate region: aT < < gT, the perturbative decay is kinematically 
suppressed due to thermal masses of quasi-particles, and a non-zero dissipation rate comes 
from their thermal widths [El [9]. As a result, the dissipation coefficient can be approxi- 
mately expressed as 



dim(r 




2/FD("^<^/2)] for > 2m^^ 



(2.12) 



+ 



for 



2m^\ < m^. 



where /fd denotes the Fermi-Dirac distribution. Aq and Ai are numerical constants, and 
in our numerical computation with a = 0.05, they are given by Aq — 0.3 and Ai ^ 2 x 10~'^ 
respectively. Note that this result is applicable to all the region in a small amplitude 
regime: A0 < rri^^ ~ gT. 



3 Reheating after inflation 

In the last two subsections, we have introduced basic ingredients to study the dynamics 
of oscillating inflaton. In this section, let us study the dynamics of oscillating inflaton and 
numerically evaluate the reheating temperature with some examples. 

3.1 Effective dissipation rate of the inflaton 

The equation of motion of the inflaton is given by 

4>+{3H + T^)^ + ml(P = 0, (3.1) 

where is the 0-dependent dissipation coefficient and H is the Hubble parameter. In 
order to study the dynamics of oscillating scalar field, it is convenient to consider quantities 
averaged over a time interval that is longer than the oscillation time scale but shorter than 
the Hubble and dissipation time scale. Taking this time average, one finds 

+ 3Hp^ = -F;%, (3.2) 

Prad + 4iJprad = rfp^, (3.3) 
3Mp2 = + p^^^^ (3.4) 



7 



where Mpi is the reduced Planck mass, prad stands for the energy density of radiation, the 
energy density of inflaton is given by := 0^/2 + m^02/2 and the effective dissipation 

rate is defined as F^^ := F^^^/^^. Here the time-averaged quantity is represented by 
Since F<^ depends on 0, F^^ is different from F,^ in generaL Let us summarize the effective 
dissipation rate F^^ relevant to our following discussion. The complete list is shown in 
Appendix [Bl 

The effective dissipation coefficient is independent of the amplitude in the the small 
amplitude regime (A0 < m^^ ~ gT) [See Eq. (1B.6P ]: 



Ff ~ dim(r) { 



X^aT 
27r2 



2 



^ J 1 - 4^ [1 - 2/fdK/2)] for 2ml < 



Ao + Ai 



aT 



for rrifp < ^m^. 



(3.5) 



For our numerical computation, we take a = 0.05, and then the numerical constants Aq 
and Ai are given by Aq — 0.3 and Ai ^ 2 x 10~^ respectively. The effective dissipation 
rate for non-perturbative production is also independent of (p and it is given by 



■^effl 



INP 



dim(r 



for X(j) ^ max 



(3.6) 



On the other hand, if the amplitude is larger than m^, then the effective dissipation 
coefficient depends on 0. In particular, the effective dissipation coefficient for the regime 
rrirf, < aT is important in the following discussion. In this case, the effective dissipation 
coefficient can be approximated as [Eq. (]B.16P 



4 ~ AT^ 
C — -Aq dim(r) — ^ for A0 > T. 

Stt-^ a 



(3.7) 



Here we have explicitly included an uncertainty in a numerical constant C that is caused 
by the extrapolation in the intermediate regime as we already mentioned, and it is taken to 
be unity: C = 1 for our numerical computation. Practically, the dissipation coefficient for 
m^jj < A0 < T is not important, since this term F^''^ oc 0^ [See Eq. ( IB. 191) ] cannot complete 
the reheating of the Universe because it decreases faster than the Hubble parameter. [See 



also Fig. d] and footnote #9 



3.2 Evolution of inflaton /plasma system and reheating 

Now we are in a position to calculate the evolution of infiaton/plasma system after in- 
flation. We numerically solve the differential Eqs. (13. 2p - (13. 4p to study the dynamics 
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Figure 1: Tlie evolution of various quantities as a function of Hubble scale H: the effective 
dissipation rate except for non-perturbative particle production T {red thin solid), one for non- 
perturbative particle production Fnp (green thin dotted), the energy density of radiation pi-ad {ma- 
genta thick solid) and inflaton (black thick dashed) normalized by an initial energy density pini- 
Top: (m^,X,(l)i) = (1 TeV, lO'^, lO^^ GeV), Middle: (m^,X,^i) = (1 TeV, 10-^ lO^^ GeV), 
Bottom: (m^. A, (/)i) = (1 TeV, 10^^ 10^^ GeV). 
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of inflaton/plasma system, using the effective dissipation rate summarized in App. IB.2I 
To make our discussion concrete, let us assume that x is an extra quark hke matter 
charged under SU(3) with the fundamental representation in the following. Then, we have 
dim(3) = 3 and T(3) = 1/2. And the gauge coupling constant is assumed to be a = 0.05 
hereafter. 

Fig. [1] shows the evolution of various quantities as a function of Hubble scale H\ the 
effective dissipation rate F^^ except for non-perturbative particle production, that for non- 
perturbative particle production F^'^l^p, the energy density of radiation p^ad and inflaton 
P(^, normalized by an initial energy density Pinii= '"^00?/2)- 

The top panel is computed with (m^. A, (pi) = (1 TeV, 10~^, 10^^ GeV). First, the radi- 
ation with high temperature (T ~ 10^ GeV) is produced via the instant preheating. The 
condition for non-perturbative production [Eq. f l2.3p ] soon saturates, since the amphtude 
scales as oc a"'^/^ where a is the scale factor of the Universe while the temperature 
scales as T oc a~^^^, and consequently the non-perturbative production shuts off. Then, 
as can be seen from the plateau of prad around if ~ 5 x 10~^ - 10~^ GeV, the tempera- 
ture of thermal plasma becomes nearly constant during the regime where the dominant 
dissipation rate is given by F^''^ ~ AT^/(a0). At this regime, the energy density of radi- 
ation behaves as prad ~ ^'fp<i>/H ~ M^^T"^ {H / (p) . Therefore, the temperature becomes 
constant since the inflaton dominates the Universe at that time. Finally, the reheating is 
completed at F^''^ ~ H. In the top panel, the reheating takes place via F^*^ ~ AT^/ (a0) at 
H 2 X 10~^ GeV, and the reheating temperature is Tr ~ 10^ GeV. 

The middle [bottom] panel is computed with (m^. A, (pi) = (1 TeV, 10~^, 10^*^ GeV) 
[{m^,X,(pi) = (1 TeV, 10^'^, 10^*^ GeV)]. The subsequent evolution is the same as the top 
panel case in the both middle and bottom panels. First, the thermal plasma is produced via 
the instant preheating, and the condition for non-perturbative production soon saturates. 
Then, the plateau region follows ii ~ 5 x 10"^ - 10"^ GeV [ii ~ 5 x 10"^ - lO^^QgV]. 
After that, since (p decreases due to the cosmic expansion, the dominant dissipation rate 
becomes F^^ ~ A^aT. In the middle panel, the reheating takes place via F^^ ~ X^aT 
at if ~ 10~^ GeV, and the reheating temperature is Tr ~ 10^ GeV. On the other hand, 
in the bottom panel, the reheating occurs via F^^ ~ A^m^ at if ~ 10~^^ GeV, and its 
temperature is given by Tr ~ 3 x 10^ GeVJ^^ 

Analytically, the reheating temperature can be roughly estimated as follows in the 

Usually, the reheating temperature is defined as the temperature at which the radiation dominated 
Universe begins and it roughly corresponds to the epoch ii ^ as In the present situation with 

thermal dissipation effect, this definition is ambiguous because of the peculiar behavior of T'^ . As seen in 
the middle panel of Fig. [l] can once become equal to H but the relation prad ^ P4> may hold thereafter 
without exponential decay of the inflaton for a while. This is because the dissipation rate decreases faster 
than the Hubble parameter during the regime: F^^ cx 0^. Therefore, the reheating temperature Tr here 
is defined as the temperature at which the inflaton energy density begins to decrease exponentially. One 
should note that, although the parameter Tr is a convenient quantity which describes a global picture 
of the early Universe, actual thermal history before the reheating would be significantly different from a 
conventional one. 
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Figure 2: The reheating temperature Tr as a function of A is shown. Left: = ITeV and 
Right: m0 = 10^ TeV. 



three cases: reheating via (i) Vf ~ XT'^/{a(j)), (ii) ~ X^aT and (iii) Ff^ ~ A^m^ : 



Ci/2 



y4odim(r) 
5'* a 



1/2 



(AMpim,^)^/' 



\ 1/2 

ylQdim(r)^a^ 
5'* 

dim(r)\^^^ .1/2 
(A^Mpim^) 



(3. 



5'* 



1/2 



(iii) 



Note that the resultant reheating temperature contains the uncertainty C from Eq. (13.71) . 
Importantly, the coupling A dependence differs among (i) - (iii) and the initial amplitude 
(pi dependence is absent even in the case (i). These behavior can be seen in Figs. [2] and [31 
In Fig. [21 reheating temperature is plotted as a function of A for m<^ = 1 TeV (left) and 
= 10^ TeV (right) with 0j = 10^* GeV. It is seen that in the small A limit, the reheating 
temperature is determined by the standard perturbative decay scenario (case (iii)). As 
Tr increases and approaches to for larger A, it begins to saturate due to the effect 
of thermal blocking. For larger A, however, thermal dissipation comes in and again Tr 
increases [case (ii) and (i)]. This figure does not depend on (pi for (pi > 10^^ GeV. In Fig.[3l 
contours of reheating temperature as a function of A and are shown. 



As mentioned in footnote \4^9\ the reheating temperature Tr here is defined as the 
temperature at which the infiaton energy density begins to decrease exponentially. The 
sharp discontinuity between two regimes [(i) and (ii)] seen in Fig. [2] is related to the 
definition of reheating temperature Tr. The reheating cannot be completed during the 
regime where the effective dissipation rate is given by F^^ oc 0^, with the definition of 
reheating that we employed. This is clearly seen in the middle panel of Fig. [H (A = 10~^): 
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Figure 3: Contour plot of reheating temperature Tr as a function of A and rrifp; Top: (pi = 
10^8 GeV and Bottom: (pi = 10^^ GeV. Inside the shaded region, the condition Xa(p > is 
violated, and this region depends on the initial amplitude (pi. At the upper left corner of bottom 
panel, one can see the region where the non-perturbative production is completely absent because 
\(pi < m^. 



12 



r crosses H two times, but at the first crossing the reheating is not completed and the 
reheating temperature is roughly determined at the second crossing. For larger A, the 
situation becomes close to the top panel of Fig.[T](A = 10~^), where reheating is completed 
soon after the first crossing. Thus the reheating temperature jumps somewhere around 
10~^ < A < 10~^ (for m0 = 1 TeV) if it is plotted as a function of A. This is the reason 
for the behavior in Fig. [2l 

Here we note that, for the purpose of estimating the reheating temperature, it is not 
necessary to impose the condition that the light degrees of freedom should be always kept 
in thermal equilibrium for every 0's oscillation. The only requirement is that the typical 
interaction time scale of plasma becomes much faster than the Hubble parameter and the 
dissipation rate of before the reheating is completed. 

4 Conclusions and Discussion 

In this paper we have investigated the issue of reheating and thermalization after inflation. 
If the inflaton is not heavy enough and its coupling to light species is not so small, the 
standard reheating scenario in which the perturbative decay of the inflaton triggers the 
reheating may not hold. In such a case, we need to carefully study the particle production, 
thermalization, and the resulting dissipation effect on the inflaton coherent oscillation. 
Actually, we found that the dissipation effect in thermal plasma plays a crucial role in 
the completion of reheating, and the reheating temperature can be much higher than the 
inflaton mass. This is consistent with the statement of Refs. |8l |9], but our setup and 
methods are more general and have broad applicability to concrete models. 

For example, in the Higgs inflation pT| \T2\ scenario (see Refs. [2T1 |22] for its real- 
ization in NMSSM), the inflaton mass around the vacuum is weak scale. Although (p'^ 
potential is dominant after inflation and the evolution of inflaton/plasma system would 
be different from our model, the final process of reheating cannot be understood without 
taking the dissipation effect into account, as studied in this paper. MSSM inflation [13] is 
another example in which the inflaton is very light. Non-perturbative particle production 
is expected at the first several oscillations [TT], but the final reheating may be caused by 
dissipative effects in thermal plasma. Similar results may hold for alchemical inflation 
scenario [H], in which the inflaton has a mass of soft SUSY breaking scale and oscillates 
around the origin after inflation. It will also be useful for reheating after a class of thermal 
inflation model [23] . It should also be noticed that the evolution of thermal plasma before 
the complete reheating in these cases can significantly differ from a conventional scenario. 
Phenomenological consequences, such as relic abundance of heavy particles created in 
high-temperature plasma and the efficiency of baryogenesis, may not be characterized by 
a single parameter Tr, but by a detailed thermal history before reheating. We will further 
generalize our results and apply to concrete models in a separate paper. 

If the inflaton oscillates around a large VEV as in the case of new inflation [2], it is 
expected that the dissipation effect would be much milder and the conventional reheating 
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scenario by the perturbative decay would be appropriate for broader parameter spaces. 
This might also be true in hybrid inflation [211 125], although it is rather non-trivial because 
the inflaton oscillates around the origin while the waterfall field oscillates around a large 
VEV. We leave these issues for a future work. 

Finally we comment on the fate of particle-like excitations of inflaton. In the case 
where thermal dissipation plays a crucial role for the reheating with Tr > m,^, the infla- 
ton particles are expected to have (nearly) thermal abundance. When the temperature 
decreases to T ~ after the completion of reheating, the inflaton freezes out from 

thermal bath. At this stage, the perturbative decay of the inflaton opens since the tem- 
perature is lower than the inflaton mass, and the decay rate is given by ~ \^m^, which is 
much larger than the Hubble rate at this epoch, H ~ lO^^m^/Mp unless A is extremely 
small. Therefore, the inflaton particles decay into light species and disappear as soon as 
they freeze out from thermal bath. 
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A Thermalization 

In this section, let us briefly review the thermalization of weakly coupled plasma for the 
case of instant preheating, following Ref. [18]. 

As discussed in Sec. 12. 2[ due to the break-down of adiabaticity of the coupled particles 
X around </> ~ 0, x particles are produced non-perturbatively in each oscillation as shown 
in Eq. (12. 4p . After the passage of origin ~ 0, they become heavier and the decay 
rate becomes larger correspondingly. Eventually, they decay into the other light degrees 
of freedom at T ^{(pit Acc))'tdcc ~ 1 before the comes back to the origin, if the decay 
rate of Xi given by P^, is sufficiently large. Assuming the typical decay rate of x as 
P^ ~ ~ 9'^aX\(j){t)\, one finds that this is the case for m</, -C O'^aXcj). Here 9 denotes 

the mixing between x and other light degrees of freedom. In this case, energy density of 
the other light degrees of freedom is given by 

Sp ~ m^nj^^^ ~ r^a-i/^A^mJ^^. ^^^^^ 
in each oscillation. The typical momentum of decay products is roughly given by 

Q-m^L^-e-'a-'/^^X'/'m'J^'/^. (A.2) 
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Therefore, the converted energy density can be expressed as 



6p ~ e^a^'/^Q^ (A.3) 
Hereafter, we assume that the mixing ^ is (9(1). In fact, the size of mixing angle is not so 



important in the following discussion. [See footnote ^11 ] 

At the first passage of ~ 0, the total energy density of light degrees of freedom is esti- 
mated as p = (5p ~ a^/^Q^. This is the so-called under occupied case since the momen- 
tum distribution / around the typical momentum Q can be evaluated as f{Q) 
It is instructive to compare the typical momentum Q with the "temperature" , defined as 
Tf ~ p^/^. This temperature can be expressed as Tf ~ a^/^Q, and hence it is smaller 
than Q. This means that the typical phase space distribution of produced particles for 
the first crossing of ~ are concentrated on the UV regime, compared to the thermal 
equilibrium distribution. As discusse d in [T8] , in this UV dominated case, the subsequent 
thermalization takes place as follows 



#10 



(i) Soft particles are radiated from the hard particles, and a new population around small 
momentum is created. They eventually fall into a thermal-like distribution below a scale 

Pmax • 

/soft(p) ~ %/P\ for p < Pinax- (A.4) 



(ii) Then, the typical scale Pmax evolves towards UV regime. The evolution of Pmax is 
dominated by the elastic scattering with the hard particles. When Pmax reaches T*, the 
distribution function becomes comparable to order 1, and then the soft sector is partially 
thermalized. This time scale can be evaluated as 

t ~ «-^/2g-i_ 5) 

At this stage, the soft sector dominates the screening effect, the number density and the 
elastic scattering. However, the energy density is still dominated by the remaining hard 
particles. 

(iii) Finally, the remaining hard particles lose their energies to the soft "thermal" sector 
by multiple splittings of daughter particles. The system thermalizes when this process is 
completed. The time scale can be estimated as 

teq ~ {a^T,)-' ^fQJTi - a~'^'''Q-\ (A.6) 

'^^^ In the following, we will only consider gauge bosons since they dominate the equilibration because 
of the induced emission factor. 
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Therefore, the equihbration time scale is determined by the time scale for the hard particle 
Q ^ T{ to lose its energy in the presence of thermal bath with temperature Tf. 



If the equilibration time scale teq, given by Eq. ( 1A.6I) . is much smaller than the oscillation 
time scale m^^, we can safely assume th at t he produced light particles have enough time 
to thermalize. This condition is given b y^^^ 

1 < a^^^/^^^A^/^. (A.8) 

B Dissipation coefficient 

In this section, we summarize the dissipation coefficient for the sake of completeness. 



B.l Definition of effective dissipation coefficient 

The equation of motion for scalar field is given by 

+ (3i/ + r<^)0 + mj^ = (B.l) 

Here T^f, is an amplitude dependent dissipation coefficient. We want to calculate averaged 
quantities with a time-interval that is longer than the oscillation period but shorter than 
the Hubble time scale and dissipation time scale. In the following, this time-average is 
represented by The energy density of (p field is defined by 



= (B.3) 

Here represents an amplitude of cj). Using the virial theorem, one can derive the evolution 
equation for the energy density: 

p^ + ?>Hp^ = -Vfp^ (B.4) 
where the effective dissipation coefficient is defined as 



Tf := (B.5) 

In general, the dissipation coefficient depends on 0, and hence the effective dissipation 
coefficient Vf has a non-trivial (j) {p^) dependence. 



"^^^ If we keep the mixing angle 9, this condition becomes 



1 < e^/'^a^'^/^y X4>/m^. (A.7) 
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B.2 List of Dissipation Coefficient 

Let us summarize the effective dissipation coefficient F? as a function of (p and m^. 



#12 



(i) A0 <C "T-^h ~ fl'^- -'■'^ ^^is case, the effective dissipation coefficient is independent of 
(p. Therefore, it is exactly the same as the dissipation coefficient r^,: 



= dim(r) < 




for m0 < 2m*^(T) 



[1 - 2/FD(m^/2)] for 2m'^{T) < m^. 



(B.6) 



Ao and Ai are numerical constants, and they are given by Aq ~ 0.3 and Ai ~ 2 x 10~^ in 
our numerical calculation with a = 0.05. Note that we neglect the hole contribution for 
simplicity |26j . 

(ii) A0 ^ m^j^ ~ gT and <^ cxT: In this case, the dissipation coefficient relevant 
to the following calculation is given by [10] 



Aq dim(r) 
Aq dim(r) 



X^aT 



for A0 ^ m^j^ 

for m^^ < A0 < T 

for T <^X(f) 



(B.7) 



where 



T(r)\^(127r)^ 



IGtt^ 



In a 



(B.8) 



where dim(r) is the dimension of x's representation r of gauge group and T(r) is the index 
of x's representation r that is defined by T{r)6"'^ = tr[t'*(r)t^(r)]. Note that the above 
dissipation coefficients are computed with two limits; small and large amplitude. Hence 
we have some ambiguities in the intermediate regime. Using these equations, one can 



"^^"^ In what follows, we consider the case where the interaction time scale in thermal plasma is much 
faster than the dissipation coefficient of inflaton: aT ^ F^. 
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compute the effective dissipation coefficient and it is given by 



27r2 



X 



sin2x 



+ 74odim(r) 



7r/2 TT 

y' sin 4y' 



sin 4x' 



+ 



02 



7r/2 
4 



27r 



7r/2 



27r 



4v 

— - 2 + 
TT TT tan y 



(B.9) 

(B.IO) 
(B.ll) 



wfiere x and y are determined case by case as follows. Aq is a numerical constant and it 
is given by Aq ~ 0.2 for our numerical computation with a — 0.05. 

(ii-i) X<p < g^T'^/m^j,: In this case, the non-perturbative production does not occur, and 
hence x, x', y' and y are given by 



m. 



arcsm ■ 



th 



X — X 



for rri^YL < 



for m^j^ > 



T 



y^y ^ 



arcsin — - for T < Xd) 
X(f) 

for T > X^, 



(B.12) 



(B.13) 



It is instructive to study the asymptotic behavior of in two cases: (a) A0 ^ T and 
(b) m^^ <^ Xcf) <^T. In the case of (a), x — x' and y — y' are given by 



m 



X — X 



y = y 



< 1, 



X(f) 
T 

< 1. 

X(f) 



Therefore, the effective dissipation coefficient can be approximated by 



4 T 



AT2 



rf ~ — Aodim(r)-^ for A0 > T. 
On the other hand, in the case of (h), x = x' and y = y' are given by 



m. 



X — X 



th 



y^y 



A0 

TT 



< 1 



(B.14) 
(B.15) 

(B.16) 

(B.17) 
(B.18) 
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Then, the dominant contribution to dissipative coefficient can be expressed as 

rf ^Aodim(r)^. (B.19) 



(ii-ii) X<p > g^T^ / m^: In this case, the non-pcrturbativc production occurs. Inside the 
region |0| < {171^0/ X)^^"^ =: </)np, the adiabaticity is broken down. Hence, we cannot use 
the WKB solutions inside this region. There are some ambiguities to evaluate the region 
where the dissipation coefficient is replaced by one caused by the instant preheating, but 
we simply evaluate the threshold value as 0np- Then, one can show that the threshold 
value 0NP is always greater than m^j^ from the inequality for non-perturbative production. 
Therefore, one finds x — and the remaining x', y' and y are given by 



0, 0, arcsin — 



W, y', y) = < 



A0, 



for 0NP > T/X 



arcsin arcsin arcsin—^ for T/A > (/)np (> Wv,/A). 
A(/) A(/. A0/ * 



(B.20) 



(in) X(f) ^ ~ gT and > T: In this case, the dissipation is dominated by the 
perturbative decay. Thus, the effective dissipation coefficient is given by 



r \2 



rf = < 



X^m^ 
Stt 



where 



Ay 4 
TT TT tan y 



y — arcsin — 



Here the coefficient c is given by 



A0 



dim(Ad) /T(r 



47r 



47r 



for A0 <^ 
for X(j) ^ 



Note that if A0 ^ m^, the non-perturbative production occurs inevitably. 



(B.21) 



(B.22) 



(B.23) 



(iv) \4> > ■ni^h ~ fl'^ ^'^'i < ''^</> < This is the missed region of our calculation. 
We simply extrapolate between (ii) and (iii) as a rough approximation. 
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